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ABSTRACT 

The ultracompact X-ray binary 4U 1820-30 is well known for its »170-d superorbital modulation in X-ray flux and spectrum, and the 
exclusiveness of bursting behavior to the low hard 'island' state. In May-June 2009, there was an exceptionally long 51-d low state. 
This state was well covered by X-ray observations and 12 bursts were detected, 9 with the high-throughput RXTE. We investigate the 
character of these X-ray bursts and find an interesting change in their photospheric expansion behavior. At the lowest inferred mass 
accretion rates, this expansion becomes very large in 4 bursts and reaches the so-called superexpansion regime. We speculate that this 
is due to the geometry of the inner accretion flow being spherical and a decreasing accretion rate: when the flow geometry nearest to 
the neutron star is spherical and the accretion rate is low, the ram pressure of the accretion disk may become too low to counteract 
that of the photospheric expansion. In effect, this may provide a novel means to probe the accretion flow. Additionally, we observe 
a peculiar effect: the well-known cessation of X-ray bursts in the high state is too quick to be consistent with a transition to stable 
helium burning. We suggest an alternative explanation, that the cessation is due to the introduction of a non-nuclear heat source in the 
neutron star ocean. 

Key words. Accretion, accretion disks - X-rays: binaries - X-rays: bursts - stars: neutron - X-rays: individual (4U 1820-30) 



1. Introduction 

Low-mass X-ray binaries with neutron stars (NSs) emit X-rays 
from two distinct locations: the inner region of the accretion 
disk, including perhaps a boundary layer and corona, and the NS. 
The former is powered by the liberation of gravitational poten- 
tial energy and the latter by thermonuclear burning in the surface 
layers consisting of accreted matter. Per nucleon the thermonu- 
clear energy is just up to a few percent of the gravitational energy 
and would be hardly distinguishable if it did not have a different 
time scale. Thermonuclear burning is mostly unstable and 'vis- 
ible' as short bursts of X-ray radiation while gravitational en- 
ergy is liberated throughout the disk on typical time scales of at 
least a few weeks. X-ray bursts are characterized by few-second 
rises and decays at least ten times as slow, and a cooling black 
body spectrum of temperature kr < 3 keV. For reviews of X- 
ray bu rsts, see iLewin et"ail d 1993b and Strohmayer & Bildsten 
(I2006h . 

Different kinds of thermonuclear X-ray bursts exist. One dis- 
tinction that one can make is the presence of photospheric ra- 
dius expansion (PRE). The luminosity produced in some bursts 
is so powerful that the Eddington limit is reached and the pho- 
tosphere is lifted fro m the NS surface dGrindlav et alj 11980; 
Ebisuzaki eTaTII 19831) . PRE can be recognized by a temporarily 
increased emission area of a cooler thermal spectrum and a bolo- 
metric flux that remains approximately constant f or the same du- 
ration . PRE happens in about 20% of the bursts (Galloway et alJ 
2008). Usually the increase in emission area is approximately 
tenfold. In 1 % of all cases, however, the increase is much larger: 
10 4 or more (or, > 10 2 in terms of photospheric radius). This is 



qual ified as 'superexpansion' din 't Zand & Weinberdl20lol see 
also lHoffman et alJI 1978trvan Paradiis et alj|1990t) . 

There is an intricate interplay between accretion and burst 
behavior (e.g.. jFuiimoto et alJ Il98lt Ivan Parad iis et all ll988: 
ICornelisse et alj|2003l) . Primarily the accretion rate determines 
the burst recurrence time: the faster the accretion and ergo fuel 
supply, the faster the bursts recur. At the highest accretion rates, 
thermonuclear burning becomes stable and bursts disappear. A 
second order effect involves the hydrogen abundance: a high 
abundance may incur a burning regime where stable hydrogen 
burning co-exists with unstable helium burning. 

The aforementioned effects are observationally obvious. It is 
to be expected, but harder to prove observationally, that there 
is also a dependence on accretion geometry. Non-isotropic ac- 
cretion on the NS surface will change the local (or specific) 
accretion rate and will influence the ignition conditions (e.g., 
Bildsten| [T998l) because those may really depend on the spe- 
cific rather than global rate. A very nice illustration of this 
is the recently discovered transient IG R J17480-2446 (e.g., 
ICavecchi et al.ll201lt iLinares et al.ll2.01 lb . 

In conclusion, accretion parameters (rate, geometry, compo- 
sition ) determine burst para meters such as the recurrence time 
(e.p JCornelissg et al. 2003), peak luminosity and duration (and, 
ergo, fluence; Gall oway et al. 2008), P RE and the amplitude of 
burst oscillations dMuno et alJl2004t IWattsI l2012h . Conversely, 
thermonuclear bursts may influence the accretion. The lumi- 
nosity of bursts is sometimes so high that it may affect the 
accreti on flow for a brief p eriod. The influence may be dra- 
matic. Kuulke rs et alJ ((2009) discuss an X-ray burst from Cen 
X-4 which appears to have turned the accretion disk into a hot 
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Fig. 1. The evolution of the 3-12 keV photon flux relative to that 
of the Crab source (black curve; RXTE-ASM data) and the 15- 
50 keV photon flux with respect to 0.18 times that of the Crab 
source (red dash-dotted curve; Swift-BAT data). Vertical lines 
indicate the bursts detected, solid for those with superexpansion 
and dashed for those without. Three bursts are indistinguishable 
because they have too short recurrence times to be visible on 
this scale (see Table [TJ. The solid horizontal lines indicate the 
days that the PCA observed 4U 1820-30 under ObsID 94090. 
The solid vertical lines in the top left corner indicate the 2<x un- 
certainties. 



state, although it appears at odds with accretion disk theory. 
Other exam ples of influence are bursts that cool the accretion 
disk corona dChen et al.ll201 ll) and bursts with superexpansion, 
such as the superburst from 4U 1820- 30 ( Strohmaver & Brown] 
120021: Ballantv ne&"Strohmaveill2004l:lKeeldl20ia and a num- 
ber of bursts from (candid ate) ultracompact X-ray binaries 
din 't Zand & Weinbergl2 010). During these bursts, episodes ex- 
ist when the flux is lower than the pre-burst accretion flux. This 
implies that accretion was stopped or covered by the optically 
thick shell. Better yet, strong factor-of-2 fluctuations were re- 
cently observed in the decay of several burst light curves that are 
very pro bably due to structural d i sturbances of the accretion disk 
surfac e din 'tZand et alj 1201 lbt iPalmerf 1201 U iDegenaar etal] 
120121) . 

Because of the influence of superexpansion on the accretion 
flow, it makes sense to study superexpansion as a function of 
accretion. Unfortunately, superexpansion is usually a rare phe- 
nomenon which makes such a study difficult. However, in 2009 
a series of observations was done on 4U 1820-30 when it was 
gradually changing its accretion rate and presumably geometry. 
Simultaneously, a series of bursts was observed with different 
levels of photospheric expansion. We here report our findings of 
these observations. 

4U 1820-30 is the sour ce in which th e X-ra y burst phe- 
nomenon was discovered by Gri ndlav et al.l dl976l) . It is a per- 
sistently bright low-mass X-ray binary (LMXB), viewed at 
an inclination angle below roug hly 50° (lAnderson etalJll997t 



iBallantvne & Strohmaverl l2004t) . It is located in the globular 
cluste r NGC 6624 at a distance of 8.4 + 0.6 kpc dValenti et al.l 
2004) and is extraordinary since it has t he shortest orbita l 
period known for any LMXB: 11.4 min (IStella et al.lll987l) . 
The orbit al period puts it in the ultracompact X-ray binary 
regime dNelson et al.lll986h . The donor stars in this group dis- 
tinguish themselves as being hydrogen-deficient dwarf stars. 
The donor in 4U 1820 -30 is probably a helium white dwarf 
dRappaport et al. 119871) . From an analysis of burst energetics, 
ICummingl (120031) conclude that the hydrogen abundance X < 
0.1. 4U 1820-30 shows a supero rbital modulation in the X- 
ray flux with a pe riod of ^170 d dPrie dhorskv & Terrell| ["l984; 
IChou & Grindlav! 120011: iFarrell et al.l 120091) . The flux resides 
mostly in a high state, except for few-day periods every » 170 
d when t he (bolometric) flux de creases by a factor of approxi- 
mately 2. IZdziarski et al.l (l2007al) studied the orbital modulation 
as a function of superorbital phase and found that the orbital 
amplitude is up to a factor of 2 lower during the superorbital low 
state (or island state), and that the orbital peak phase shifts. This 
suggests a different geometry of X-ray emitting and scattering 
structures on the accretion disk between both states. The nature 
of the w 170 d periodicity is speculat ed to be the result of the 
influ ence of a third massive body, see |Chou & Grindlav (2001) 
and lZdziarski e t al. (2007b), which modulates the orbital eccen- 
tricity of the ultracompact system and, thus, the mass accretion 
rate. 

It is only during the low states that the source shows X-ray 
bursts. This is attributed to the difference in accretion rate: in the 
high state thermonuclear burning is stable and not distinguish- 
able from the accretion radiation while in the low state it is un - 
stable and visible as X-ray bursts (e.g JChou & G rindlav 2001b . 
4U 1820-30 is the only confirmed ultracompact hydrogen- 
deficie nt X-ray binary to have exhibited carbon-fueled s uper- 
bursts (Strohmaver & Brown 20021 liiT't Zand et al.ll2.Ql lal see 
also IBallantvne & Strohmaverll2004l) . 

The subject of our study is the bursting behavior during the 
longest low state ever observed in 4U 1820-30, lasting 51 days 
in May- June 2009. 

2. Observations 

Figure Q] shows two light curves of the long low state. The 2-12 
keV RXTE-ASM d ata (black; for a description of the ASM, see 
iLevine et al.| [T996) show the typical low/high state behavior. At 
MJD 54938 the flux sharply decreases from 0.30 Crab units to 
0. 10 Crab five days later. This endures until MJD 54994 (5 1 days 
later) when there is a sharp increase lasting 4 days and reaching 
0.40 Crab. During the low state, the flux shows a slow bowl- 
shaped change with a minimum halfway through of 0.06 Crab. 
The bowl is a little skewed, as though there is a slow increase 
over a time scale longer than the low state. 

The 15-50 keV BAT data (red curve, as obtained from the 
Swift/BAT Hard X-ray T ransient Monitor web sitJt for an in- 
strument description, see iBarfhelmv et all 12005) show more or 
less opposite behavior, with slowly monotonically increasing 
fluxes during the ASM low state, from 0.05 Crab to 0.16 Crab at 
the end of the ASM low state. 

Twelve bursts were detected during the long low state with 
three instruments. They are listed in Table [T] Three were de- 
tected with BAT, of which one triggered a slew that resulted in a 
detection of the next burst with XRT and BAT. Nine bursts were 
detected with the Proportional Counter Array (PCA). Three pairs 
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Table 1. Specifications of all detected bursts detected during the long low state (# 1-12) and 2 bursts that have the smallest hard 
color value in the RXTE archive (# A-B; see Fig . . Numbers between parentheses refer to the uncertainty in the last digits. 



Id. Instr. 


Date Time MJD RXTE OB SID 


PCUs 


SC 


HC 


Fluence 


Decay Pers.flux 4 EED 5 












(10~ 7 ergcnT 2 ) 


time 3 


(10~ 9 erg 


(s) 




CUT) 








meas. 


extrap. 


(s) 


cm - s ) 




1 PCA 


Apr 27, 2009 19:42:36 54948.82125 94090-01-01-02 


2 


1.92 


1.08 


3 2(2) 


4.1 


3.3(2) 


4.01(3) 


0.70 


2 PCA 


Apr 29, 2009 16:52:04 54950.70283 94090-01-01-05 


1.2 


1.96 


1.14 


3.3(2) 


4.3 


3.4(1) 


3.93(3) 


0.90 


3 PCA 


May 05, 2009 18:35:37 54956.77474 94090-01-02-03 


2.4 


2.05 


1.19 


3.4(2) 


4.5 


3.2(1) 


4.11(3) 


0.90 


4 PCA&BAT 


May 07, 2009 17:45:34 54958.73999 94090-01-02-02 


2 


2.04 


1.20 


3.4(2) 


4.3 


3.4(2) 


4.04(3) 


1.30 


5 PCA 


May 27, 2009 07:42:57 54978.32149 94090-01-04-00 


2 


2.08 


1.24 


3.5(2) 


5.0 


3.7(2) 


4.75(4) 


1.68 


6 PCA 


May 27, 2009 11:52:39 54978.49490 94090-01-04-01 


2 


2.03 


1.25 


3.6(2) 


5.2 


3.6(2) 


4.81(4) 


1.55 


7 PCA 


May 30, 2009 04:29:42 54981.18730 94090-01-05-00 


0,2,4 


2.12 


1.24 


3.6(1) 


4.6 


3.5(1) 


5.00(2) 


1.40 


8 BAT 


May 31, 2009 21:38:16 54982.90157 


















9 BAT 


lun 06, 2009 09:28: 19 54988.39465 


















10 XRT&BAT Jun06, 2009 12:43:37 54988.53029 


















11 PCA 


Junl2, 2009 12:49:14 54994.53420 94090-02-01-00 


2 


1.92 


1.12 


3.2(2) 


4.2 


3.1(2) 


6.73(4) 


0.59 


12 PCA 


lun 12, 2009 14:42:45 54994.61303 94090-02-01-00 


1,2 


1.95 


1.09 


3.2(1) 


3.6 


3.5(1) 


6.73(4) 


0.41 


A PCA 


Sep 29, 2004 10:29:03 53277.43685 90027-01-03-05 


0,2,3,4 1.74 0.97 


3.5(2) 


4.6 


3.2(1) 


4.28(2) 


0.61 


B PCA 


Mar06, 2011 18:32:11 55626.77236 96090-01-01-020 


1,2 


1.74 0.92 


2.4(2) 


2.9 


2.7(2) 


4.90(2) 


0.10 



1 This is the fluence as measured under the exponential decay fitted to the time profile of the bolometric flux after the peak, plus the peak flux 
times the time interval between the peak flux and the start of the precursor; 2 This is the fluence as measured under the exponential decay of the 
bolometric flux, backwards extrapolated to the start of the precursor; 3 The e-folding decay time as determined from the bolometric flux time 
history from the peak value on; 4 2-60 keV flux from fit with comptonized model to all data from the relevant observation. Superexpansion 
equivalent duration (see text) 



of bursts had short wait times of 4.2 hr (bursts 5 and 6), 3.3 hr 
(bursts 9 and 10) and 1.9 hr (bursts 1 1 and 12). Interestingly, the 
last burst was detected on the same day as the low state ended 
with the sharp increase in ASM flux (Fig. [TJ). The moment of 
change was not covered by PCA observations. 

The nine bursts detected with the high-throughput PCA on 
RXTE provide much better statistical quality than the others. 
Therefore, we focus on these. During the low state, 4U 1820- 
30 was observed 27 times for a total of 134 ks with the RXTE 
PCA (ObsID 9 4090). Figure [ TJ show s the coverage. 

The PCA dJahoda et alJ 20061) consists of 5 proportional 
counter units (PCUs) that are sensitive between 2 and 60 keV 
and combine to give an effective area of 6000 cm 2 at 6 keV. 
The spectral resolution is about 18% full width at half maxi- 
mum at 6 keV and the time resolution of the data products for 
4U 1820-30 is typically 125 /is. It is rare that all PCUs are active 
at the same time, particularly at late times in the mission, but 
the central PCU number 2 (counting from 0) is almost always 
on. Table [TJindicates which PCUs were on per burst. A wide va- 
riety of data products can be requested from the PCA. For our 
analyses, we used standard-1 (0.125 s resolution and 1 channel 
for 2-60 keV photons) and standard-2 data (16 s resolution and 
128 channels) and event mode data E_125us_64M_Q_ls (125 /is 
resolution and 64 energy channels). 

Fig.[2]shows the light curves of all PCA bursts as drawn from 
the standard-1 data of all operational PCUs. It is immediately ob- 
vious that for bursts 4, 5, 6 and 7 the photon flux drops below the 
pre-burst level, resulting in the appearance of a precursor imply- 
ing superexpansion. The photosphere expands while emitting at 
the Eddington limit. The increase in emission area is compen- 
sated by such a strong decrease in radiation temperature that the 
thermal spectrum shifts out of the PCA bandpass. Additionally, 
the accretion radiation is cut short by the expanding shell, either 
sweeping the inner accretion disk or covering the X-ray emitting 
part of the disk. We note that the XRT burst (number 10) is the 
first after burst 7, after six days, with sufficient statistical quality 
to see superexpansion if it were present, but none was detected. 



The 4 superexpansion bursts are annotated by solid vertical lines 
in Fig. QJ (three pairs of bursts are indistinguishable in the plot 
because they occur within a few hrs). 

We include in our analysis 2 other bursts from the total of 17 
detected by the PCA in its 16 year mission (for a discussion of 
14 of these, see lKeekll2012h . because they are the bursts with the 
least amount of expansion. The top two light curves in Fig.[2]are 
of these 2 bursts. 



3. Analysis 

3. 1 . Time-resolved burst spectral analysis 

We divided each of the 1 1 PCA bursts in a number of time bins 
for which we extracted spectra. The time bins were chosen such 
as to be able to follow spectral changes promptly. A spectrum 
of the non-burst emission was extracted from standard-2 data in 
a time interval between 160 and 16 s prior to each burst. The 
non-source background (cosmic diffuse emission or particle- 
induced) was calculated for each extracted spectrum through 
pcabackest version 3.6 and taken into account in subsequent 
analyses. The non-burst source emission was modeled by a sim- 
ple spectral shape (thermal bremsstrahlung) that adhered to the 
data. Response matrices were calculated for each observation us- 
ing pcarsp version 11.7.1. Allowance was made for a system- 
atic error of 1 % per channel. 

The burst spectra were modeled by a Planck function for 
black body radiation leaving free the temperature and emission 
area. Allowance was made for the non-burst emission to vary by 
including the spectral model as obtained for the pre-burst inter- 
val, with all parameters fixed, and leaving free its normalization. 
From the fitted Planck function, bolometric fluxes were calcu- 
lated. For PRE phases, it is to be expected that the shape of the 
non-burst emission will change. We do not model this and keep 
it fixed. Therefore, the black body parameters are expected to be 
biased for these phases and we do not draw conclusions from 
them. We did not include low-energy absorption in our models, 
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Fig. 2. Measured photon rate time histories for the 9 PCA bursts, 
from bottom to top chronologically, supplemented with the 2 
PCA bursts that are lowest in the color-color diagram (Fig. @). 
The subtracted background level is measured from a 150 s time 
interval prior to the burst and the curves have been normalized 
to the peak rate. For guidance, the zero level has been indicated 
by a dotted line. 



because the column density (7Vh ~ 1 x 10 21 cm 2 ; Kuulker s et alj 
l2003t ICostantini et al.ll2012l) is too low to be noticeable with the 
PCA. Fig. [3] shows an example of the fit results, for a burst that 
was measured with 3 PCUs (number 7, see Table[TJ. 

The e-folding decay time r was measured by fitting to the 
time history of the bolometric flux, beyond the peak, an expo- 
nential decay function. The 'precursor time' was determined as 
the time between the peak of the flux and the start of the pre- 
cursor. The 'gap time' was estimated by investigating the light 
curve and estimating the time between the end of the precursor 
decay and the rise of the main burst phase. The burst fluence 
was determined by combining the integral of the fitted exponen- 
tial decay function with the assumed fluence before the peak flux 




0.10 1.00 10.00 

Time [sec] 

Fig. 3. Time-resolved spectroscopy of burst number 7. The top 
panel shows the measured full-bandpass photon count rates for 
the PCA Xenon layers (solid line; units of 10 4 s" 1 ) and propane 
layers (dashed line; units of 3000 s _1 ). The propane layer has a 
somewhat softer res ponse and th erefore better tracks the lower 
energy photons (e.g., Keek 2012). 



which is simply determined by the peak flux times the precursor 
time. It should be noted that this measures only the nuclear en- 
ergy that is transformed into radiation. The nuclear energy may 
also transform into kinetic and gravitational energy of expelled 
material. A rough estimate of the energy lost to expelled material 
may be obtained by back-extrapolating the exponential decay to 
the time of the start of the precursor. We find that this amounts 
to between 10 and 40% of the (presumed) total nuclear energy. 
The values are largest for the four superexpansion bursts, mostly 
because the precursor times are longest there. The results of the 
measurements are in Tableland Fig. [5] 

There is little variation in burst fluence and duration. The val- 
ues vary between 3.2 and 3.6xl0~ 7 erg cirT 2 , a 12% range with 
a relative lcr uncertainty of about 5%, and is consistent with be- 
ing constant (x 1 = 11.4 with 8 d.o.f. for a fit with a constant). 
The weighted average is (3.39 + 0.05) x 10~ 7 erg crrT 2 . For a 
8.4 + 0.6 kpc distance, the fluence translates to a total energy of 
£burst = (2.9 + 0.4) x 10 39 erg. If this is due to a pure layer of he- 
lium, with a 30- energy production of E nuc = 5.84 x 10 17 erg g _I , 
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Fig. 4. Color-color diagram of all PCU2 measurements of 4U 
1820-30 betwee n March 31, 1999, a nd March 31, 201 1. Colors 
are defined as in lHoman et al.l (120101) : soft color is the count rate 
in PCU2 channels 10-17 to 6-9 (corresponding to ranges 4.1- 
7.4 and 2.5-4.1 keV) and hard color is channels 24-40 to 18-23 
(17.1-9.9 and 7.4-9.9 keV). Each point is for 128 s of data. Burst 
data have been excluded. The stars refer to the colors of the 128- 
s data stretches immediately before the 4 superexpansion bursts, 
the open circles to those of the other 5 bursts in the low state, 
squares to all other ordinary bursts and the large cross to the 
superburst. One ordinary burst is not included here. The burst 
recurrence time is 2.5 times larger at the group of circles with 
coordinates (1.95,1.09) than at the group of stars with coordi- 
nates (2.05,1.25). 

the ignition column depth is j/j gn « Zsburst/47r/? 2 £ nuc = (4.0 ± 
0.7) x 10 8 g cirT 2 for a canonical NS radius of 10 km. The upper 
limit on the variation in fiuence (10%) is, unfortunately, not ac- 
curate enoug h to constrain the hydrogen abundance in a manner 
suggested by Cummin iH (l2003h . 

3.2. Spectral modeling of accretion-powered radiation 

To put the data i n perspective, we s how in Figure [4] the com- 
monly used (e.g.,|Homan et al] |2010l) color-color diagram of all 
PCA data of 4U 1820-30 obtained after March 22, 1999. We 
exclude earlier data (including one burst) because early gain 
changes of the detectors preclude straightforward inference of 
colors. The data show two branches, commonly referred to as 
the banana state (lower branch ) and island state (upper branch) 
( Hasi neer & van der Klisl fl989i. All low states are equivalent to 
the island state. Also shown in Fig. [4] are the data points for the 
17 PCA bursts detected in this period. These are taken from the 
128-s intervals immediately prior to the bursts. Obviously, all 
bursts occur in the island state. Most notably, the superexpan- 
sion bursts occur in the extreme upper part of the island state. 

It is difficult to find an unambiguous physical model for the 
spectra of the accretion radiation, or 'persistent spectrum' as 
we shall call it. Low-resolution LMXB 2-30 keV spectra are 
most often modeled with a combination of (disk) black body and 
comptonization components, representative of the thermal radi- 
ation from the accretion disk, NS and/or boundary layer, and the 
comptonization in a corona in the neighborhood of that disk. The 
geometry of the corona is uncertain. We tested this model against 



the spectra of the complete observations for each burst, employ- 
ing PCA standard product spectra. In particular, we employed 
a simple black body component a nd a Comptonize d component 
according to the prescription by Titarchuk (1994). In XSPEC 
they are identified by model names 'bbody' and 'comptt'. 
We find that the evidence for the black body component is very 
uncertain in these data and that including it in the model de- 
grades the accuracy with which comptonization parameters can 
be measured, due to strong coupling between the two compo- 
nents (both are very broad continua without unambiguous dis- 
tinctive significantly detected spectral features). Therefore, we 
choose to ignore the black body component. This implies some- 
what worse fits, but the xl values remain fairly close to accept- 
able - between 1 and 3.9. The spectra are well represented by 
the comptonization component within a level of 3% per channel. 
We here present the fitted parameters, with errors that have been 
multiplied by a factor of yjxl which would be equivalent with 
a fit with xl — 1 ■ The resulting 2-60 keV fluxes are provided in 
Table Q] and Fig. [5] plasma temperatures and optical depths (for 
a disk geometry) are provided in Fig. [5] We note that the en- 
ergy of the seed photons is always far below the lower bandpass 
threshold and, thus, not constrained by the data. The spectral re- 
sults show that the plasma temperature and optical depth varies 
significantly by a factor of 3 for 2-60 keV fluxes that are compa- 
rable. 

For a 8.4 kpc distance, the 2-60 keV fluxes translate to a 2- 
60 keV luminosity range of 3.3 to 5.7 x 10 37 erg cm _2 s _1 . For a 
canonical 10 km/1 .4 M NS this is equivalent to a mass accretion 
rate of 2.8 to 4.8 xl0~ 9 M Q yr or a specific mass accretion rate 
of 1.4to 2.4xl0 4 g s-'cirT 2 . 

3.3. Correlation study 

We introduce a new empirical parameter to measure the amount 
and duration of expansion: the equivalent expansion duration 
(EED). The basis for this is the background-subtracted normal- 
ized light curve as shown for all bursts in Fig. [2] The EED is 
similarly defined as the equivalent width of absorption lines in 
spectra: it is the loss of photons due to the expansion, expressed 
in seconds that it would take to emit those photons at the level 
of the peak flux. The EED is a function of the bandpass and the 
effective area function, but this is of no consequence here be- 
cause we use the same instrument in its full bandpass through- 
out. When superexpansion occurs, the light curve drops to nega- 
tive values. This will push the EED to even higher values. 

We plot various burst parameters against other burst pa- 
rameters as well as accretion parameters in Fig. [5] The only 
strong one-to-one correlation exist between EED and the shape 
of the pre-burst spectrum as measured through the comptonizing 
plasma parameters. 

4. Discussion 

The data can be summarized as follows: 

• the fiuence and decay rate of all bursts are equal within 6%, 
except for burst B. The helium content is presumably close 
to 100%, see below. This implies that the difference in photo- 
spheric expansion is unlikely to be the result of a difference 
in nuclear energy generated or ignition depth; 

• the amount of expansion is strongly correlated with the po- 
sition of 4U 1820-30 on the color-color diagram in the sense 
that the 4 bursts with the strongest expansion (superexpan- 
sion) are at one extreme end of the island state and those with 
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Fig. 5. Diagrams of burst parameters versus persistent emission 
parameters (left column, except bottom plot) and burst param- 
eters versus burst parameters (right column). The plotted flu- 
ence is the measured fluence (see Table[TJ. The dashed blue data 
points refer to bursts A and B (top 2 light curves in Fig. [2]). The 4 
superexpansion bursts are plotted with thick red lines. For addi- 
tional information a diagram is shown between plasma tempera- 
ture and optical depth (bottom left plot). 



the smallest expansion at the opposite end. Equivalently, the 
EED is strongly correlated with the comptonization parame- 
ters; 

This strongly suggests that the accretion environment deter- 
mines the amount of expansion exhibited by thermonuclear X- 
ray bursts. The question is: in what way? Is it related to the rate 
of accretion, the geometry of the accretion, or both? 

All changes in the persistent radiation are thought to be 
driven by changes in accretion rate. By how much does the 
accretion rate change? It is well known that fluxes do not di- 
rectly translate to accretion rates. It appears to be true only on 
order-of-magnitude scales. One can think of various reasons why 



there is no 1-to-l correspondence between flux and accretion 
rate: the measured persistent flux is not bolometric, the flux may 
have changing anisotropy factors and the radiation efficiency 
of transforming gravitational to radiation energy may change. 
Fortunately, our data provide an independent means of measur- 
ing M changes. The three recurrence times measured for X-ray 
bursts from 4U 1820-30 (see Sect. in combination with the 
burst fluences, indicate that the accretion rate at burst 12 (with- 
out superexpansion) is 2.5 ±0.2 times larger than at burst 6 (with 
super expansion; provided that all nuclear burning is unstable 
and that the nuclear and gravitational energy liberated per ac- 
creted nucleon is constant; i.e., the burst alpha parameter, see 
iLewin et~aDll993L is constant). For comparison, the 2-60 keV 
persistent flux ratio is only 1 .4+0. 1 over these bursts. Note that 
this ratio difference supports t he pure hel i um ac cretion hypothe- 
sis: the pure helium model in lCummind (l2003h predicts that M 
should be 1 .5 to 2 times higher than the value derived from the 
X-ray flux to have agreement with burst recurrence times and 
energies in 4U 1820-30. 

Changes in spectral shape (as measured via for instance color 
in Fig. |4]or plasma temperature in Fig. [5]) are generally attributed 
to changes in the geometry of the inner accretion flow, proba- 
bly as a result of changes in accretion rate. The low/hard and 
high/soft states are very common in variable LMXBs, disregard- 
ing the nature of the a ccretor (NS or black hole; see reviews 
by Ivan der Klisl l2006t iRemillard &"McClintock1l2006l respec- 
tively). These distinct states are usually accompanied by distinct 
Fourier spectra. 

PCA timing and spectral studie s of 4U 1820-30 have 
been the subject o f previous papers. IZhang etaO d!998l) and 
Blose r etaD (|2000) report on twin kHz quasi-periodic oscilla- 
tions (QPOs) and find that they, from the lower left banana 
state to the island state, graduall y move in tandem to frequen- 
cies twice as small (see also lAltamirano et alj 2005 ). In the 
banana state , the f requencies saturate (but see Mendez 120021) . 
Blos er et alj (|2000) modeled the PCA spectra by a combination 
of a black body and comptonization component and find that 
in the island state the plasma temperature increases and the op- 
tical depth decreases. We observe the same trend further into 
the island state (Fig. |5). This behavior is consistent with an in- 
terpr etation that the QPO s originate in the inner accretion disk 
(e.g., Ivan der Klil[2006) and that the inner edge of the accre- 
tion disk moves outward (recedes) from the lower left banana to 
the upper island state and that the cleared region is filled by a 
non-thermal more spherical advection-dominated (i.e., with less 
radiati ve efficiency) accretion flow. This is supported by the find- 
ings of Zdziarski et al. (2007a) mentioned in Sect.Q] 

We note that we have searched for QPOs in the new data, 
but were able to find only a single, broad QPO around 540 Hz 
in the combined power spectra of the three observations with a 
soft color lower than 2.0 (i.e. roughly middle of the island state). 
A comparison with the p ower spectra of 4U 1820-33 analyzed 
in lAltamirano et al.1 d2005l) suggests that this QPO is the upper 
kHz QPO. We note that our data are statistically less constrain- 
ing than those of Altamirano et al. (2005) due to a smaller num- 
ber of active PCUs. The other five observations, with soft colors 
higher than 2.0, extend to more extreme parts of the island state, 
in which the inner accretion disk edge possibly moves out even 
further. No high frequency QPO is found in these observations, 
although indications for a broad (Q 0.6) feature around 280 Hz 
are seen. 

Expansion into the circumstellar medium will become easier 
when the ram pressure of the accretion flow becomes smaller. 
This may happen if the accretion rate decreases or if the geom- 
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etry changes. The latter may occur, for example, if the accretion 
is channeled through a disk instead of spherically. A large part of 
the solid angle as seen from the NS will then open up and allow 
free expansion. However, all bursts that we observe occur during 
the island state, when the inner flow is thought to be spherical all 
the time and the inner edge of the accretion disk is thought to be 
at some distance from the NS. If the frequency is Keplerian in 
nature, a factor of 2 decrease in frequency (see above), implies 
a factor of 2 2/3 = 1.6 increase in radius from banana to island 
state. 

The change in expansion factor is, therefore, most likely to 
be directly related to the accretion rate. In the following we at- 
tempt to find quantitative support for this. The ram pressure of 
free f alling spherical accretion can be expressed as dFrank et al.l 
119921 page 135) 



(1) 



with p the density of the infalling matter and vg the free-fall ve- 
locity 



v s = {2GM/R) 1 ' 2 

The density can be derived from the continuity equation 

4nR 2 pv = M 

so that Eqn.[T]becomes 

(2GM) 1/2 . s/ , 
— M R~ 5/2 

M /T 5/2 dyne cirT 2 



(2) 



(3) 



An 

= 1.5 xlO 12 



(4) 



for a 1 .4 M G NS, with M in g s 1 and R in cm. 
The ram pressure of a shell leaving the NS is 



Aram.shell - P Expand 



yo 
i 



R 



expand 



(5) 



where yo is the initial column thickness of the shell, / the geo- 
metrical thickness and u e xpand the expansion velocity of the shell. 
The shell will be able to leave the NS surface (R — /?ns) if 



ram, shell 



>P ra , 



M < 6.7 x 10" 



or if 



Vo 2 R 

i "expand ^NS 



5/2 



(6) 



Once the shell has left the NS, it depends on how the shell thick- 
ness evolves whether its ram pressure will remain above that of 
the accretion flow. For a constant shell thickness it does, but for 
a shell thickness that grows stronger than R 5 it does not. 

What values can be applied for the parameters in Eq. [6]? 
We use the initial rate of increase in the black body radius (see 
Fig. [3]) as a rough estimate of the expansion velocity: i> exp and ~ 
10 3 km s _1 . The shell column thickne ss t/n may be derived 
from the superexpansion duration (see lin 'tZand & Weinberg! 
120101) . It takes about 1 s for the shell to become transparent 
(at a column thickness of 1 g cm" 2 ) due to dilution that goes 
as R 2 . During that time the shell has expanded to a radius 
of ^expand x( ~ 10 3 km. Back extrapolating the dilution im- 
plies that yo ~ 10 4 g cm" 2 . This is 10~ 4 times the ignition of 
?/ign ~ 4 x 10 8 g cm" 2 . For the geometrical layer thickness one 
may use the pre ssure scale height of a helium atmosphere which 
is ~ 33 cm (e.g. JPiro & Bildstenl2007h . For# NS = 10 6 cm, Eq.|6] 
predicts that the shell will be able to lift off if M ■£> 3 x 10" 5 



M e yr" 1 . This is a very loose constraint which is fullfilled under 
any practical circumstance, 10 4 times larger than measured for 
4U 1820-30. This may be due to the relative arbitrariness of the 
chosen parameter values. For instance, I may increase substan- 
tially during the expansion and u e xpand ma y evolve. If Z is -100 
m instead of 33 cm and t> e xpand 3 times smaller, the constraint 
in Eq. [6] would be similar to the observed M. Furthermore, we 
have used a very simple approach (isotropic spherical free-fall 
accretion, isotropic shell expansion, constant shell velocity and 
thickness, ignoring radiation pressure, and ignoring GR correc- 
tions). The conclusion is that a quantitative verification of our 
proposed scenario of the superexpansion behavior in 4U 1820- 
30 needs detailed modeling which is outside the scope of the 
present paper. 

The constant ignition depth would imply that the expelled 
shell thickness is similar in all bursts, but that the amount of 
distance traveled by the shell differs. Usually, it does not travel 
far, perhaps less than 100 km, before falling back. Only for the 
superexpansion bursts does it travel far enough to move out of 
the X-ray band. This is supported by the fact that there is al- 
ways evidence for a shell, even if there is no superexpansion: all 
bursts from 4U 1820-30 show a drop in flux within 1 s from the 
burst onset. This idea is supported by the slower decay of the ini- 
tial spike for bursts without superexpansion. The shell is moving 
slower and does not reach high altitudes (i.e., in Fig.[2]the spikes 
drop the fastest in the superexpansion case). 

We note that the superburst in 4U 1820-30 also has super- 
expansion while it is on the opposite side of the island state in 
the color-co lor diagram (Fig. 0J. A new analysis of this event 
(lKeekll20T2l) shows that the energy contained in the precursor 
is larger than those in ordinary bursts from 4U 1820-30, which 
may be explained by a shock induced by the carbon detonation 
and reaching the photosphere. Therefore, the driving of the ex- 
pansion occurs differently from that in ordinary bursts from 4U 
1820-30 and the shell ram pressure is likely to be different. 



5. A note on the lack of X-ray bursts in the banana 
state 

As an addendum, we address a peculiar phenomenon in our mea- 
surements: X-ray bursting activity seems to react more rapidly to 
changes in the accretion-induced radiation than expected from 
crustal heating. The first burst after the flux arrives in the low 
state is detected within 5 days. The last burst before the persis- 
tent flux returns back to the high state is detected only two days 
before. No bursts were detected in PCA measurements within 
the next 8 days with a total exposure of 10.4 ks while the re- 
currence time between the two last bursts was 6.8 ks s, nor with 
Swift observations. We made a similar but more significant mea- 
surement in 201 1, when we followed up a PCA burst detection 
with a 1.2 d long continuous Chandra observation starting 1.3 d 
after the burst. The jump to the high state occurred in the 1.3 d 
time frame between both observations. No burst was detected in 
RXTE data and Chandra data within 2.5 d after the burst with a 
total exposure time of 1 .4 d. In the same time frame there were 
0.3 d of exposure time with Swift-BAT, without any X-ray burst 
detection. We note that a similarly sud den stop of bursting a ctiv- 
ity was observed in IGR J17473-2721 (IChenevez et al.ll201 lh . 

The data suggest that bursting activity stops as suddenly as 
the low state does. This is p eculiar. In the common picture of 
4U 1820-30 if X < 0.03 (ICummingl 120031 the bursting ac- 
tivity is driven by the temperature of the crust. The crust is 
heated through pycnonuclear reactions and electron-capture re- 
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actions (e.g jGupta et al.ll2007l:lHaensel & Zdunikll2008l) whose 
rates are a direct function of the accretion rate. However, 
the thermal time scale o f the crust is of order months (e.g., 
Brown & Cumming 2009), so there should be a delay between 
a change in the accretion rate and one in the bursting activity. 
Therefore, the prompt lack of X-ray bursts in the high state of 
4U 1820-30, and vice versa, is not explained by stabilization of 
helium burning by crustal heating. 

An alternative explanation might be that in the banana state 
the inner edge of the accretion disk is probably closer to the NS, 
see Sect. [4] Depending on the inclination angle and thickness of 
the inner disk and/or boundary layer, this is perhaps sufficient 
to block our view of the NS and make X-ray bursts invisible. 
A test of this hypothesis would be the detection of partly ob- 
scured X-ray bursts with a few-hours recurrence times. This has 
never been observed in 4U 1820-30, despite abundant coverage 
throughout the history of X-ray astronomy. Furthermore, X-ray 
bursts have been observed in banana states o f other sources (e.g., 
lvanderKliset~ai1ll990t iMuno et alJl2000h . Therefore, this ex- 
pla nation is unlikel y. 

IChenevez et al.l (1201 lb proposed another explanation for the 
sudden stop of bursting activity in IGR J17473-2721: the occur- 
rence of a superburst during a data gap. However, given the long 
recurrence time of superbursts (of the order of 1 yr for persistent 
sources and longer for transients such as IGR J17473-2721), this 
appears to us an unlikely scenario. 

A more promising explanation seems that during the high 
state the primary heating of the flash layers does not occur by 
the crust, but by a heat source at a shallower depth. A thermal 
time scale of ^ 1 d applies to column depths of ;£ 10 12 g cm -2 
dBrown & Cum ming 2009) which is in the NS ocean and out- 
side the crust. The heat source is most likely not of nuclear ori- 
gin. For example, it cannot be CNO burning, because there is 
no significant abundance of hydrogen in the accreted material. 
Inogamov & Sunvaev (2010) suggest that for accretion rates in 
excess of a few percent of the Eddington limit, rotational heat- 
ing may stabilize helium burning. In this model, the braking of 
a hypersonic azimuthal flow in the accretion disk by the NS at- 
mosphere results in an additional heat source. Perhaps the ac- 
cretion flow geometry in the banana state of this source allows 
for more efficient transfer of angular momentum from the disk 
to the NS envelope, boosting this heat source and quenching 
the bursts. Quantitative assessment of this hypothesis requires 
2-dimensional multi-zone simulations of nuclear burning, pri- 
marily as a function of NS spin, accretion composition and flow 
geometry. Unfortunately, the NS spin has not been measured yet 
for 4U 1820-30. 
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